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Abstract Vortical flows in shallow water interact with long surface waves by 
virtue of the nonlinear terms of the fluid equations. Analytical formulae are 
derived that quantify the spontaneous generation of such waves by unsteady 
vorticity as well as the scattering of surface waves by vorticity. In a first Born 
approximation the radiated surface elevation is linearly related to the Fourier 
transform of the vorticity. The "dislocated" wavefronts that are analogous to the 
Aharonov-Bohm effect are obtained as a special case. 
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Nonlinear effects in one-dimensional surface wave propagation over a shal- 
low fluid have been widely studied in the past two decades in association with 
solitonsjl]]. What happens in two dimensions is, by comparison, far less understood! 
The effect of vorticity and of submerged bodies 0] on free surface motion has 
also been studied, again in one dimension, and the vortex-surface interaction is 
the subject of much current research 

In a different vein, in recent years it has become apparent that the nonlinear 
interaction between sound and vorticity at low Mach numbers can be profitably 
understood in terms of concepts borrowed from classical field theory of sources 
and waves |]]. This has allowed for an understanding of vortex dynamics in a 
slightly compressible fluid and has suggested the use of ultrasound as a probe 
of vorticity in both ordered and disordered flows[H]- Recent experimental results 
have confirmed the soundness of this proposal|]]. On the other hand, it is well 
known [|ToH that the nonlinear propagation of sound obeys similar equations as the 
nonlinear propagation of surface waves in shallow water, and it is the purpose 
of this paper to apply the concepts that have been of use to understand the 
interaction of sound with vorticity to the understanding of the interaction of two- 
dimensional surface waves with ordered or disorderd vorticity in shallow water. 

Consider an incompressible fluid of undisturbed uniform depth h moving with 
velocity v in a uniform gravitational field g. We shall refer to it as shallow water 
although of course it need not be water. The free surface will be described by 

z = C(?,y,t) 

where ( = h + £(x,y,t), and ( is the deviation of the surface away from the 
horizontal, whose typical length scales will be supposed to be much longer than 
h, thus allowing for the neglect of surface tension. In shallow water vertical 
variations of v± are neglected and the governing equations are 

^ + v ± -(c^±) = o 

r\ — * 

-^ + (v L -V ± )v ± = ~gV ± ( (1) 

where the subindex _L means horizontal, or x — y, components and will be omitted 
from here on. Equations ([!]) are supposed to be evaluated at z — (, and viscosity 
has been neglected. The boundary condition at the bottom is then v z (z = 0) = 
0. These equations are similar to, but not identical with, the equations for a 
compressible, adiabatic bulk flow and they can be combined to yield 

- ghV\ = hV ■ ((v-V)v) - |(V • (vQ). (2) 
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This is a wave equation for the surface waves with a source term due to the 
nonlinear couplings. One can think then of at least two situations of interest 
when there is a bounded (two dimensional) vortical flow with typical velocities 
very small compared to yfgh: one is the spontaneous generation of surface waves 
by this flow, in analogy with the aeroacoustic generation of sound by vortical 
flows. The other is the scattering of surface waves by the flow. 

The reason for the surface wave generation by a bounded flow is that, although 
surface deformations will not be the dominant effect at the source, far away from it 
they will decay, because of their wavelike nature, like (distance) -1 / 2 , while the two 
dimensional flow associated with the source, being incompressible, will decay like 
(distance) -1 . A formal solution to Eqn. (pj) can be written as a convolution of the 
right hand side with the Green function for the two dimensional wave equation. If 
surface elevations are neglected at the source, an acceptable approximation, the 
source looks almost like the Lighthill source termfjH]]. The similarity is, however, 
not an identity because in the present case Vj_ • v±_ ^ 0. 

In the spontaneous generation of surface waves by a flow of bounded (in space) 
vorticity the second term in the right-hand-side can safely be neglected because 
surface deformations are supposed to be negligible at the source and one obtains 
the following expression for the far field surface waves spontaneously radiated by 
such a flow: 

( = A(cJam)* vg 

cr 

where u is a perfectly two dimensional incompressible flow and Q = V A u is 
its (two-dimensional) vorticity. This is in very close analogy to the formulae for 
aeroacoustic sound derived by Powell and Howe [|T^ 1 which are best manipulated 



in Fourier space 13 , where 



G(x,v) = —H+' l/l1 ' 1 



57T ° V C J 

is a Hankel function with outgoing wave boundary condition. 

The case of scattering of a plane wave can be studied with ease in the case 
that wave frequencies v are high compared to any frequencies Q associated with 
the target flow, and fluid velocities v s associated with the surface wave are small 
in comparison with those associated with the vortical flow, u, which in turn are 
supposed to be small compared with c = \fgh. To this end we write 

v = v s + u (3) 

where u is a perfectly two dimensional incompressible vortical flow with a per- 
fectly flat surface and v s (v a <C u) is what is needed for v to satisfy Eqns. ([!]). 
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Moreover, the differential equation (0) can be turned into the integral equation 

C = Cinc + G*S (4) 

where d nc is the incident plane surface wave and G * s is a convolution of the 
Green function G for the wave equation with the source 

s = /iv.((tr.v)tO-^(v.(trc))- (5) 

Substitution of the decomposition (|) into the source © and neglecting the 
horizontal components of the vorticity, a valid approximation in shallow water, 
as well as terms quadratic in v s , leads to a source that is the sum of three com- 
ponents: 

s = Si + s 2 + s 3 

where 

S\ = hV ■ (Co x v 



S2 = h(\7 2 (u-v s )-^- t (u-V)() 
s 3 = W-((m-V)m). 

The source S3 corresponds to the spontaneous generation of surface waves 
discussed above and it will not be considered further in the scattering context 
since the waves so generated are of much lower frequency than the incoming and 
scattered waves. For weak surface waves the scattered amplitude will be much 
weaker than the incident one and Equation (ffl) can be solved in a first Born 
approximation, that is replacing v s by the incident plane wave value 

Wine = v n cos(fc • x - u t), 
for which velocity at the surface and surface elevation are related by 

Using the fact that, in the far field 

_ _ x dG 
c at 

and integrating by parts it is tedious but straightforward to show that if the 
scattered wave is written as 

Cscatt = Cl + C2 



4 



where Ca — G * s a , a — 1, 2, then 



The second term on the right is of order ~ Quuvq/c 2 and can be neglected with 
respect to ( 2 , which is of order ~ k^uvo, in the case under consideration of 
incident frequencies much higher than typical frequencies associated with the 
vortical target: Q <C u Q . 

We are thus finally led to the relation 

cos# 

Cscatt = — f n T\ G * (V • (O; X -y inc )) (6) 

g{cosu — 1) 

where 9 is the scattering angle and Q the vorticity. Using the far-field expression 
forG, 

v ' ' 4tt \ 2ttvx / c 4 J 



this expression becomes 



which is the promised relation between scattered surface wave amplitude at po- 
sition x and frequency v in terms of to, the Fourier transform of vorticity. The 
incident wave has (velocity) amplitude v and frequency z/ ; q is the momentum 
transfer. Note that the apparent divergence of this expression for small scattering 
angles is an artifact of taking an incident plane wave which is, technically, infinite. 
The angular dependence will break down for wavelengths A at angles such that 

sin 6 ~ X/L 

where L is a length of the order of the spatial extent of the incident wave. Formula 
(0) suggests a possible non intrusive way to probe vortical flows in shallow water 
using surface wave scattering. 

Berry and collaborators 0] have remarked that a plane wave of wave vector k 
propagating along the x direction incident on a stationary point vortex will give 
rise to surface deformations that can be locally (as opposed to globally) described 
by a multivalued phase: 

£ ^ ^i(kr cos 0+vt+aO) 
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corresponding to wavefronts exhibiting a dislocation at the position of the vortex. 
Here (r, 9) are polar coordinates centered at the vortex and a = uT/c 2 , where T 
is the vortex circulation. For the dislocation to be noticeable a must be of order 
one. It is of interest to see how this is related to the results reported above, which 
are not immediately applicable since the scattering of surface waves by vortical 
flows was worked out for high frequency waves, for which this last condition is 
violated. Indeed, at a distance ~ A from a point vortex the frequency of the 
vortical flow is Q ~ T/X 2 which is ~ v when a ~ 1. 

In the presence of a steady divergenceless background flow U giving rise to 
surface deformations Co, small velocity v and surface deviations ( will obey the 
equations, easily derived from (jl|), 

(^ + ^' V ) C + W ^ = ° 

(JL + tj.Vjt + gC = (8) 

where <fi is the velocity potential for the velocity, v = V0, which exists locally 
outside the vortex core. Quadratically small terms have been neglected and 
position is supposed to be sufficiently far away from the vortex core that gradients 
are dominated by wavevectors. Under the additional assumption that U 2 <C c 2 , 
equations @ yield 

Note now that the velocity field outside a vortex is the gradient of a multivalued 
scalar potential which is proportional to the polar angle centered at the vortex: 

U = 

Looking now for time harmonic solutions in the form 

d> = che iut e^ 




we have that with 



obeys the equation 



2txc 2 



(v 2 + %)<p = 0, 
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giving rise to the dislocated wavefronts of Berry et. al.JT4 . 

In addition to the dislocated incident wave there is a scattered wave. In the 
case of the Aharonov-Bohm effect |]T4| the latter can be calculated using vanishing 
boundary conditions at the "vortex" core. The scattering mechanism in the fluid 
case is different, its origin being not in impenetrable boundary conditions but in 
the nonlinear interaction terms that were neglected above. The computation of 
this effect offers an interesting challenge which is, however, outside the scope of 
the present paper. 

To conclude, we have studied the nonlinear interaction between surface waves 
(for a two dimensional surface) and vorticity in shallow water. Analytical for- 
mulae for the generation and scattering of such waves by any (two-dimensional) 
vortical flow have been derived under the following assumptions: a) Time scales 
are such that viscosity can be neglected; b) Length scales are such that surface 
tension can be neglected; c) For scattering, the frequency of the incident wave is 
high by comparison with the (inverse of) the time scale of the vortical flow, and 
the particle velocity associated with the wave is supposed to be small compared 
with the velocity of the vortical flow, which in turn is supposed to be small com- 
pared to y/gh, the phase velocity of the waves. Dislocated wavefronts analogous 
to the Aharonov-Bohm effect have been obtained as a special case. The question 
naturally arises as to what happens when the fluid is not shallow. The point of 
view presented in this paper, namely that of studying the nonlinear interaction 
between surface motions and vorticity by way of succesive approximations using 
ideas borrowed from classical field theory, may be of use in studying this problem. 

This work was supported by DTI Grant E2854-9244, FONDECYT Grant 
1265-91 and CEC Contract CIl*-CT91-0947. 
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